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ABSTRACT 

We evaluate the effect of inhomogeneity energy on the expansion rate of the universe. Our method is 
to expand to Newtonian order in potential and velocity but to take into account fully nonlinear density 
inhomogeneities. To linear order in density, kinetic and gravitational potential energy contribute to 
the total energy of the universe with the same scaling with expansion factor as spatial curvature. In 
the strongly nonlinear regime, growth saturates, and the net effect of the inhomogeneity energy on 
the expansion rate remains negligible at all times. In particular, inhomogeneity contributions never 
mimic the effects of dark energy or induce an accelerated expansion. 

Subject headings: cosmological parameters — cosmology : theory — large-scale structure of universe 



1. INTRODUCTION 

Recent observations of type-la suoernov ae llRies 
and the cosmic microwave background ( B ennett et alJ 
2003) in tandem suggest that the cosmological expansion 
is accelerating. Understanding the source of this accel- 
erated expansion is one of the greatest current unsolved 
problems in cosmology (Peebles & Ratrall2003|) . Accel- 
eration seems to render inadequate a universe consisting 
entirely of matter, and appears to require an additional, 
unknown type of energy (dark energy, perhaps realized as 
a cosmological constant). An alternative to dark energy 
is that acceleration arises from a known component of 
the universe whose effects on the cosmic expansion have 
not been fully examined. One possibility currently being 
examined is that inhomogeneities in a matter dominated 
universe, on either sub-horizon (Rasanen 2005; Notari 



2005H or super- horizo n scales (jKolb et al.ll2004 l2005albl 



Barausse et al.l 120051) . may influence the expansion rate 



at late times. The central idea is that the energy induced 
by inhomogeneities leads to additional source terms in 
the Friedmann equations, with effects on the dynamics 
that leave no need for a separate dark energy compo- 
nent. In their entirety, these proposals present conflict- 
ing claims and a general state of much confusion: does 
the inhomogeneity energy produce an accelerated expan - 
sion, acting in effect as dark e nergy (Kol b et alJI2005bD . 
or does it behave as curvature l)GeshnizlaTiiTt*Tnl2Ti05() ? 
Is the magnitude of the effect small, large, or e ven di- 
vergent, on either large sc ales lIKolb et al.ll2005bj) . or on 
small scales at late times (lNotarill2005|) ? 

Part of the confusion arises from the fully relativistic 
perturbation theory formulation of many of these calcu- 
lations. Although this is undeniably a valid approach, 
the number of terms in a perturbation theory calcula- 
tion can be large and can mask the underlying physics. 
In this Letter, taking advantage of phenomenological re- 
sults that have been derived from a combination of quasi- 
linear perturbation theory, nonlinear theory, and numer- 
ical simulations, we compute the potential and kinetic 
inhomogeneity energies within the horizon to Newtonian 
order in potential and velocity for fully nonlinear density 
contrasts. We find these energies to be small at present, 
and their projected values remain small, even far into the 



future. Section 2 considers the effect of inhomogeneities, 
for weak gravity and slow motions but for arbitrary den- 
sity perturbations, characterized in terms of the density 
power spectrum. Section 3 presents the results for the 
kinetic and potential energies in both the linear and the 
fully nonlinear regimes, as a function of the cosmological 
expansion factor. Finally, section 4 discusses the impli- 
cations of these results. 

2. EFFECTS OF INHOMOGENEITIES 

The purpose of our work is to investigate whether in- 
homogeneity energy can mimic the effects of dark energy 
for a universe containing only matter. To this end, we 
work in an Q m = 1 Einstein-de Sitter universe, with 
no curvature or cosmological constant, and compute the 
effects of inhomogeneities on the cosmic expansion rate. 
The dynamics of cosmological expansion are governed by 
the Friedmann equations, 
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Any mass or energy density that makes up a significant 
fraction of the total can influence the evolution of the cos- 
mological scale factor a(t). A contribution to the energy 
density of the universe with equation of state pi = w P i 
has pi oc a~ 3 ( 1+w \ or pi/p m oc a~ 3w ; in particular, a 
component with p oc a~ 2 behaves as w = — | or cur- 
vature, and a component with constant p behaves as a 
cosmological constant or dark energy. 

We introduce t he effects of inhomo geneities following 
the formulation of Seliak & Hui ( 199*3) ■ In the conformal 
Newtonian gauge, with metric 

ds 2 = a 2 (r)[-(l + 2^)dr 2 + (1 - 2<j>)dx 2 }, (2) 

the time-time Einstein equation (G° ) yields 

3(-) 2 (l - 2V0 + (2 + 60)1 V 2 0+ 1(V0) 2 

= ^Gp{l + 5){l + v 2 ), (3) 

Where ~ ip from the space-space components of G M „. 
(Our numerical factors are slightly different from those of 
Seljak & Hui; they make little difference in the results.) 
The source on the right-hand side includes a density per- 
turbation S = Sp/p in the material rest frame, with the 
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transformation to the cosmological frame expanded to 
leading order for small v 2 . Ignoring 4>V 2 (f>, (V0) 2 , and 
v 2 , the homogeneous part of this equation reproduces 
the usual Friedmann equation. The inhomogeneous part 
reveals that (f> obeys the Poisson equation with source 
4irGpa 2 5. The volume average of the entire equation 
then leads to 

& = Y Gf5{1 ~ 5W + 2Kh (4) 

where W and K are the Newtonian potential and kinetic 
energy per unit mass, 

W=\{(l + S)<j>), K=\((l+5)v 2 ). (5) 

These expressions are correct to first order in <j> and v 2 , 
but neither an assumption nor an approximation in 5. 
We assume that ( V 2 </>) = 0; in all other places the Pois- 
son equation is adequate to determine 4>. 

The Newtonian potential and kinetic energies thus can 
influence cosmological expansion. We can compute both 
W and also K completely and exactly from knowledge 
only of the density power spectrum. The potential is re- 
lated to the density inhomogeneity by the Poisson equa- 
tion, V 2 </> = AirGpa 2 S, an expression which holds even 
for nonlinear inhomogeneities. From this, we obtain 



W 



1 



4nGpa 2 



d 3 k P(k) 



dk 



an 



2 ~"~ r ~ J (2tt) 3 k 2 J k 

expression correct in both linear and nonlinear 
regimes if P(k) is the appropriate linear or nonlinear 
power spectrum. The last equality defines the dimen- 
sionless spectral density A 2 v (k). 

In linear perturbation theory, valid for small inhomo- 
geneities, the density contrast grows as S = So(x)D(t), 
where in a mat ter dominated universe D(t) tx a(t) oc i 2 / 3 
l|Peehleslll98(H) . The kinetic energy follows from the lin- 
earized equation of continuity, 6 + V • v/a = l|Peeblesl 
1980). 

1 2 f d 3 fc P{k) 

(the usual factor /(fi) ~ ft 6 = 1 for Vt m = 1). The ki- 
netic energy scales with a(t) as d 2 D 2 , while the potential 
energy scales as pa 2 D 2 ; and so both W and K grow as 
D 2 la oc a(t), or pjj = p(W + K) cx a -2 . As was noted 
by iGeshniziani et aLl l)2005ft for super-horizon inhomo- 
geneities, in perturbation theory inhomogeneity energy 
has the same effect on the expansion rate as spatial cur- 
vature. We note that -Kii n /|Wn n | = H 2 /4irGp — §, a 
fixed ratio in the linear regime. The full kinetic energy 
in principle involves higher order correlation functions 
and is not a simple integral over the power spectrum. 
Nonetheless, the full kinetic energy can be obtained sim- 
ply from th e potential en ergy through the c osmic energy 
equation of llrvinel i|1961f) and iLavzerl l) 1963ft . 

'' 2ir '- ±)W, (8) 
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with initial conditions set in the linear regime, Kn n = 
§ | Wiin | ■ Equations JSJ) and JBJ provide us with expres- 
sions sufficient to calculate nonperturbative contribu- 
tions to the expansion rate for both the gravitational 
potential perturbation and kinetic energy components. 
The results of these calculations are given in the next 
section. 



3. RESULTS 

Equations |JB} and JSJ determine the inhomogeneity en- 
ergy of the universe as a function of epoch, which we 
characterize by the expansion factor a/a . For the pri- 
mordial power spe ctrum, we use t he CD M power spec- 
trum as given by iBardeen et alJ 1)1986). with spectral 
index n = 1, Q m = 1 , and COBE norma lized ampli- 
tude S H = 1.9 x 10" 5 l|Bunn fc White! 11997ft . To obtain 
the nonline ar power spectrum we use the lin ear-nonlinear 
mapping of iPeacock fc Doddsl lfl99i Il99^1 . The results 
of these calculations are shown in Figures 1 and 2. 
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Fig. 1. — Spectral density of gravitational potential energy 
A^r(fc) [the integrand of eq. JgJ], evaluated at the present, plotted 
as a function of wavenumber k. The dashed line shows A^, in 
linear perturbation theory; the solid line shows the fully nonlinear 
form. 

Figure shows the dimcnsionless spectral density of 
gravitational potential energy A 2 v (k) defined in eq. @), 
evaluated at the present, plotted as a function of 
wavenumber k. The dashed curve shows the density in 
linear perturbation theory, and the solid curve shows its 
fully nonlinear form. 
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Fig. 2. — Fractional contributions of gravitational potential 
energy W (long-dashed line) and kinetic energy K (solid line) to 
the total energy density of the universe, plotted as a function of 
past and future expansion factor for an f2 m = 1 universe. The 
short-dashed line is the sum of contributions from inhomogeneities. 
The dotted lines show results from linear perturbation theory. 
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Figure shows the contributions of potential energy 
and kinetic energy to the energy density of the universe, 
for past and future expansion factors in an fi m = 1 uni- 
verse. At early times, perturbation theory gives an ac- 
curate result, but at a/a rj 0.05 (redshift z w 20) the 
behavior starts to change, for an interval growing faster 
than a 1 with the fastest growth as a 1 ' 2 , and then satu- 
rating and growing significantly more slowly, eventually 
as In a. 

4. DISCUSSION 

In this Letter we have evaluated the size and the time 
evolution of the contribution of inhomogeneities to the 
expansion dynamics of a matter-dominated universe, in- 
cluding the effects of fully nonlinear density inhomo- 
geneities. When density fluctuations are in the linear 
regime, the ratio of the inhomogeneity contribution to 
the matter density grows linearly with expansion factor, 
as does curvature in an open universe, making only a very 
small contribution to the expansion rate. As density fluc- 
tuations begin to go nonlinear, the inhomogeneity energy 
grows at a slightly faster rate, at most as a 12 oc a~ 3w , 
or w — —0.4. This by itself, even if the dominant en- 
ergy component, would be only temporarily and only 
very slightly accelerating, with deceleration parameter 
<7o = 5(1 + 3iu) = —0.1. Since, at this time, the total 
fraction of inhomogeneity energy is flu w 10~ 5 <C 1, this 
has a negligible effect on cosmological expansion dynam- 
ics. 

As the universe further evolves, so that the main con- 
tributions to W and K come from deeply nonlinear 
scales, we compute the potential energy from integration 
of the nonlinear power spectrum, and obtain kinetic en- 
ergy from the cosmic energy equation (eq. In a scale- 
invariant model with power spectrum P ~ fc" as k — > 0, 
the kinetic and potential energies K and W scale with the 
expansion factor as a 

(i_„)/ (3 +n) ||D av i s fc PeeblesllT977l) 
(logarithmically in a as n — > 1), with ratio 

K - 4 

W\ ~ ? + n' ( ' 

Numerical simulations show that this continues to hold 
for the CDM spectrum with effective index n = 
d log P/d log k at an appropri ate scale, the basis of 
the l inear-nonlinear mapping llPea,cock Pod d si ITflfll 
1996). For the CDM spectrum, with n — > 1 on large 
scales, this means that growth stops, and the ratio tends 
to the virial value K/\W\ — > \ at late times. We note 
that aside from the integration of the Layze r-Irvine equa- 
tion, m any of these results were obtained bv lSeliak &: Huil 

Our results show that the contributions of the potential 
and kinetic energies of inhomogeneities has never been 
strong enough to dominate the expansion dynamics of 
the universe. For a universe with Q m = 1 today, nor- 
malized to the large scale fluctuations in the microwave 
background, the net effect of inhomogeneities today is 
that of a slightly open universe, with £1^ « 10~ 4 in cur- 
vature. The maximum contribution comes from scales 
of order 1 Mpc, falling off rapidly for smaller and larger 
k, as illustrated in Figure ^ The behavior on asymp- 
totically small scales (k 3> 10 6 /iMpc -1 ) depends on an 



extrapola tion that ignores su c h det ails as star forma- 
tion, but Fukueita fc Peeblesl 1)20041) estimate that the 
net contribution of dissipative gravitational settling from 
baryon-dominated parts of galaxies, including main se- 
quence stars and substellar objects, white dwarfs, neu- 
tron stars, stellar mass black holes, and galactic nuclei, 
is in total 10~ 4 9 of the critical energy density. 

The suggestion that nonlinear effects for large inhomo- 
geneities may mimic the effect of dark energy is not the 
case for the fully nonlinear theory. It is true that higher 
order terms in perturbation theory grow faster; the gen- 
eral n-th order term grows as D n (t) . There indeed comes 
a scale in space or an evolution in time where the behav- 
ior of higher order terms appears to diverge. Neverthe- 
less, the fully nonlinear result is well behaved. It is only 
the perturbation expansion that breaks down, and the 
actual energy saturates and grows more and more slowly 
at late times. As illustrated in Figure [3 the nonlinear 
potential and kinetic energies remain small compared to 
the total matter density at all times, even an expansion 
factor of 10 3 into the future. Inhomogeneity effects do 
not substantially affect the expansion rate at any epoch. 
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Fig. 3. — The expected fluctuation in the potential energy per 
unit mass ((AW) 2 ) 1 ^ 2 evaluated at the present as a function of in- 
frared cutoff fc m i n for n = 0.95, n = 1, and n = 1.05 (solid lines, top 
to bottom). Dashed lines are analytic approximations that asymp- 
totically become k~ °- 025 ) (log fc) 1 / 2 , or constant, respectively. The 
dotted line shows the result for a rolling spectral index that has 
n = 0.95 on the horizon today but approaches n = 1 as fc — + 0, as 
predicted by most models of slow-roll inflation. The mean value 
{ W) = 3.1 X 10~ 5 is shown as the horizontal dashed line. 

It has been pointed out that although the average in- 
homogeneity energy is small, its variance has a logarith- 
mically divergent contribution fr om the variance of the 
potential on super- horizon scales l|Kolb et alJl2005a|) . 

((AW) 2 ) = ±Jd 3 x d 3 x' ~ p ( cj>(x) cf>(x') ) 

= (27rGpa 2 ) 2 J ^L^W 2 (kR), (10) 

windowed over the horizon volume (for calculational con- 
venience we use a Gaussian rolloff rather than a sharp 
radial edge). For n — > 1 as k — > 0, this is indeed loga- 
rithmically dependent on the low-fc cutoff (and if n < 1 
the divergence is worse), but the rest of the integral is fi- 
nite for the CDM spectrum. The fluctuation in potential 
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energy, ((AW) 2 ) 1/2 , is shown in Figured] as a function 
of the infrared cutoff fc min . The integral is dominated by 
the smallest values of fc, where perturbations are deep in 
the linear regime. For n = 1 the result is very accurately 
AW = 1.45 x 1(T 5 llnk^Rnl 1 / 2 . (We note that for 
n — > 1 the units of fc m in are unimportant.) The fluctua- 
tion is comparable to the mean ( W) = 3.1 x 10 -5 when 
the cutoff is near the scale of the horizon k — Hq/c, 
and does not become of order 1 until k m i n ~ 10~ 170 
(for n = 0.95), or fc min ~ 10~ 10 (for n — > 1), or ever 
(for n > 1). While such an exponentially vast range of 
scales may not be beyond the range of possibility in an 
inflationary universe, it requires a fearless extrapolation 
well beyond what is known directly from observation. 
The fluctuation AW is dominated by contributions from 
modes that are deep in the linear perturbation regime, 
and scales with expansion factor as AW oc pa 2 D, con- 
stant in time. This contribution to the energy will ap- 
pear dynamically in the Friedmann equation as another 
matter component. Furthermore, in the presence of a 
true dark energy component, any effects on cosmologi- 
cal expansion arising from inhomogeneities quickly be- 
comes unimportant once dark energy becomes dominant 
(|Seliak fc Huilll99l . 

The fact that fluctuations in the potential diverge re- 
mains troublesome. It has been recognized for some time 
that potential fluctuations in the standard model with 
n — > 1 are logarithmically divergent, but since for most 
purposes the value of the potential is unimportant, this 
has not been perceived as a significant problem. The ef- 



fect of potential on the expansion dynamics is real, but 
the weak logarithmic divergence and the fact that it is 
a feedback of a gravitational energy on gravitational dy- 
namics may lead one to hope that this divergence is al- 
leviated in a renormalized quantum theory of gravity. 

We have found that, to leading order in <p and v 2 
but with fully nonlinear density fluctuations, inhomo- 
geneities on sub-horizon scales have only a minimal ef- 
fect on the cosmological expansion dynamics, even far 
into the future, and in particular never result in an ac- 
celerated expansion. Other authors have also shown 
that recent attempts to explain an accelerated expan- 
sion through s uper-horizon perturbations face signifi- 
cant difficulties ( iFlanae an 2005: iGeshniziani et alJl2005t 
iHirata k Seliakll2005|) . The possibility that a known 
component of the universe may be responsible for the 
accelerated expansion remains intriguing. However, we 
conclude that sub-horizon perturbations are not a viable 
candidate for explaining the accelerated expansion of the 
universe. 



We thank Dan Chung for helpful discussions on the 
general topic of cosmological inhomogeneities and Uros 
Seljak for pointing out the Seljak & Hui reference and 
for helpful comments. E.R.S. acknowledges support from 
the University of Florida's Alumni Fellowship program. 
This research has made use of NASA's Astrophysics Data 
System. 



REFERENCES 



Barausse, E., Matarrese, S., & Riotto, A. 2005, astro-ph/0501152 
Bardeen, J. M., Bond, J. R., Kaiser, N., & Szalay, A. S. 1986, ApJ, 
304, 15 

Bennett, C. L., et al. 2003, ApJS, 148, 1 

Bunn, E. F., & White, M. 1997, ApJ, 480, 6 

Davis, M., & Peebles, P. J. E. 1977, Ap JS, 34, 425 

Flanagan, E. E. 2005, hep-th/0503202 

Fukugita, M., & Peebles, f. J. K. 2004, ApJ, 616, 643 

Geshnizjani, G., C hung. D. J. H., & Afshordi, N. 2005, 

astro-ph/0503553 
Hirata, (J. M., & Seljak, U. 2005, astro-ph/0503582 
Irvine, W. M. 1961, Ph.D. thesis, Harvard University 
K olb, E. W.. Mata rrese. S., Notari, A., & Riotto, A. 2004, 

astro-ph/0410541 
Kolb, H. W., Matarrese, S., Notari, A., & Riotto, A. 2005, 

Phys. Rev. D, 71, 023524 
K olb, E W., M atarrese, S., Notari, A., & Riotto, A. 2005, 

hep-th/0503117 



Layzer, D. 1963, ApJ, 138, 174 
Notari, A. 2005, astro-ph/0503715] 

Peacock, J. A., & Dodds, S. J. 1994, MNRAS, 267, 1020 

Peacock, J. A., & Dodds, S. J. 1996, MNRAS, 280, L19 

Peebles, P. J. E. 1980, The large-scale structure of the universe 

(Princeton: Princeton University Press) 
Peebles, P. J. E., & Ratra, B. 2003, Reviews of Modern Physics, 

75, 559 

Rasane n, S. 2005, ArXiv Astrophysics e-prints, 

|arXiv:as~ ro-ph/0504005 
Riess, A. G. 2000, PASP, 112, 1284 

Seljak, U., & Hui, L. 1996, ASP Conf. Ser. 88: Clusters, Lensing, 
and the Future of the Universe, 88, 267 



